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pH ' We define Letac-Wesolowski-Matsumoto-Yor (LWMY) functions as decreasing functions from 

^0 , (0,oo) onto (0, cxj) with the following property: there exist independent, positive random 

variables X and Y such that the variables f{X + Y) and f{X) — f{X + Y) are indepen- 
dent. We prove that, under additional assumptions, there are essentially four such functions. 
The first one is f{x) — 1/x. In this case, referred to in the literature as the Matsumoto- 
Yor property, the law of X is generalized inverse Gaussian while Y is gamma distributed. 
In the three other cases, the associated densities are provided. As a consequence, we obtain 
a new relation of convolution involving gamma distributions and Kummer distributions of 
type 2. 
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Many papers have been devoted to generalized inverse Gaussian (GIG) distributions since 
their definition by Good [5] (see, e.g., [1, 8, 15, 16]). 

The GIG distribution with parameters yUCM, a,fe>0 is the probability measure 



H' ^T^/ ,^,, : fbY X^' ^ -(a^x-^+b^x)/2. 



Cp_. GIG(M,a,6)(dx-)=(^-j 2j^ {ahf +" "^/'l(o,oo)(2;)dx, (1.1) 

where K^ is the classical McDonald special function. 

(1) We stress the close links between GIG, gamma distributions and the function 
/o(x) = l/x(x>0). 

(a) The family of GIG distributions is invariant under /q: we can easily deduce 
from (1.1) that the image of GIG(/i,a,6) by /g is GIG(--/i,6, a). 



This is an electronic reprint of the original article published by the ISI/BS in Bernoulli, 
2012, Vol. 18, No. 1, 119-136. This reprint differs from the original in pagination and 
typographic detail. 



1350-7265 © 2012 ISI/BS 



A.E. Koudou and P. Vallois 



(b) BarndorfF-Nielsen and Halgreen [1] proved that 



GIG{-fi,a,b)*j(fi,—\=GlG{fi,a,b), fi,a,b>0, (1.2) 

where 7(^,6^/2) (da;) = 2i'^r{tj.) ^^~^ "^^P ~ T^^(o,oo) (x) dx. 

Therefore, if X ^ GIG(— A, a, a) and Y r^ j(X, a^/2) are independent random variables, 
then 

X^^fo{X + Y). (1.3) 

Letac and Seshadri [8] proved that (1.3) characterizes GIG distributions of the type 
GIG(-A, a,a). 

(c) Almost sure realizations of (1.2) have been given by Bhattacharya and Waymire [3] 
in the case fj,= 2, Vallois [16] for any fi> hy means of a family of transient diffusions 
and Vallois [15], theorem on page 446, in terms of random walks. 

(2) The so-called Matsumoto-Yor property is the following: let X and Y be two 
independent random variables such that 

X^GIG(--/i,a,&), r^7(M,&V2), (/^,a,6>0). (1.4) 

Then, 

U:=^^ = .UX + Y), y:=l-_l_ = /o(X)-/o(X + y) (1.5) 

are independent and 

U^GlG{~^i,h,a), F-7(Ai,aV2). (1.6) 

The case a = b was proven by Matsumoto and Yor [11] and a nice interpretation of 
this property via Brownian motion was given by Matsumoto and Yor [12]. The case 
^ = — i of the Matsumoto-Yor property can be retrieved from an independence property 
established by Barndorff-Nielsen and Koudou [2] (see [7]). 

Letac and Wesolowski [9] proved that the Matsumoto-Yor property holds for any fj,, 
a,b>0 and characterizes the GIG distributions. More precisely, consider two independent 
and non-Dirac positive random variables X and Y such that U and V defined by (1.5) 
are independent. There then exist /x,a,6> such that (1.4) holds. 

The starting point of this paper is to study the link between the function /o : a; 1— > 1/a; 
and the GIG distributions in the Matsumoto-Yor property. 

Obviously, the Matsumoto-Yor property can be re-expressed as follows: the image 
of the probability measure (on R^) GIG(— /i,a,6) (g) j{fj,,b'^/2) by the transformation 
Tfg : {x,y) i-> {fo{x + y),fo{x) — foi^ + y)) is the probability measure GIG(— pt,6,a) ^ 
7(/x,a^/2). This formulation of the Matsumoto-Yor property, joined with the Letac and 
Wesolowski result, leads us to determine the triplets {fix , IJ-y , f) such that: 

(a) fixTfJ-Y are probability measures on (0,00); 
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(b) / : (0, oo) — !• (0, cx)) is bijective and decreasing; 

(c) if X and Y are independent random variables such that X ^ fj,x and Y ^ ^y , 
then the random variables U = f{X + Y) and V = f{X) — f{X + Y) are independent. 

Unfortunately, we have not been able to solve this question without restriction. Our 
method can be applied provided that / is smooth and ^x and fxy have smooth density 
functions (see Theorem 3.1 for details). After long and sometimes tedious calculations, 
we prove (see Theorem 2.2) that there are only four classes, J"i, . . . , Jv4, of functions / 
such that Tf keeps the independence property. Then, for any / £ J^i, 1 < i < 4, we have 
been able to give the corresponding distributions of X and Y and the related laws 
of U and V (for J^2,-7^3 and J^4, see Theorems 2.4, 2.14 and Remark 2.5). The first 
class, J-"i = {a/x; a > 0}, corresponds to the known case f = fo- This case, as mentioned 
in Remark 3.3, allows us recover, under stronger assumptions, the result of Letac and 
Wesolowski that the only possible distributions for X and Y are GIG and gamma, re- 
spectively. The proof of Lctac and Wesolowski is completely different from ours since 
the authors make use of Laplace transforms and a characterization of the GIG laws as 
the distribution of a continued fraction with gamma entries. We have not been able to 
develop a proof as elegant as theirs because, with / = /o, we have algebraic proper- 
ties (e.g., continued fractions), while these properties are lost if we start with a general 
function /. 

It is worth pointing out that one interesting feature of our analysis is an original char- 
acterization of the families of distributions {/3a{a,b,c);a,b,a > 0,c £ M} and the Kum- 
mer distributions {K'-'^\a, b, c);a, c> 0,b eR} (see (2.14) and (2.29), respectively). The 
Kummer distributions appear as the laws of some random continued fractions (see [10], 
page 3393, mentioning a work by Dyson [4] in the setting of random matrices). 

As by-products of our study, we obtain new relations for convolution. For simplicity, 
we only detail the case of Kummer distributions of type 2: 

K^^'>{a,b,c)*j{b,c)^K^^'>{a + b,-b,c). (1.7) 

Obviously, this relation is similar to (1.2). 

Inspired by the result of Letac and Wesolowski [9] and Theorem 2.6, we can ask (for the 
purposes of future research) whether a characterization of Kummer distributions could 
be obtained via an "algebraic" method. 

As recalled in the above item (c), there are various almost sure realizations of (1.2) 
and of the convolution coming from the Matsumoto-Yor property. One interesting open 
question derived from our study would be to determine a random variable Z with distri- 
bution K^'^' {a + b, —b, c) which can be decomposed as the sum of two explicit independent 
random variables X and Y such that X ^ K^^^^ (a -I- b, —b, c) and Y ^ 7(6, c). 

The paper is organized as follows. We state our main results in Section 2. In Section 3 
we give a key differential equation involving / and the log densities of the independent 
random variables X and Y such that f{X + Y) and f{X) — f{X + Y) are independent 
(see Theorem 3.1). Based on this equation, we prove (see Theorem 3.9) that there are 
only four classes of such functions /. The theorems stated in Section 2 are proved in 
Section 4; however, one technical proof has been postponed to the Appendix. 
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2. Main results 

Definition 2.1. Let f: (0,oo) — > (0,oo) be a decreasing and bijective function. 

(1) We consider the transformation associated with f 

T^:(0,«))2^(0,oof, 

ix,y)^if{x + ij),fix)-fix + y)). (2.8) 

The transformation Tf is one-to-one and if f^^ is the inverse of f , then 

{Tf)-^=Tf-.. (2.9) 

(2) Let X and Y be two independent and positive random variables. Let us define 

(c/, V) = r/(x, Y) - {f{x + y), f{x) - f{x + y)). (2.10) 

/ is said to be an LWMY function with respect to (X, Y) if the random variables U 
and V are independent, f is said to be an LWMY function if it is an LWMY function 
with respect to some random vector {X,Y). 

One aim of this paper is to characterize LWMY functions. Let us introduce 

/i(2;) = -^, a;>0, (2.11) 



1' 



gi(x) = /ri(x)=ln(^i±^). 



x>0 (2.12) 



and, for S > 0, 



fUx)^\og[ ' : , ^ 1, x->0. (2.13) 



Theorem 2.2. Let f : (0, cxd) — > (0, oo) be decreasing and bijective. Under some additional 
assumptions (see Theorem 3.1, (3.7) and (3.8)), f is an LWMY function if and only if 
either fix) = f , fix) = ^fii^x), fix) = ^gii/Sx) or fix) = i/|(/3x) for some a,P, 
5>Q. 

Remark 2.3. (1) The four classes of LWMY functions are Fi = {a/x;a > 0}, J-2 = 
{i/i(/3x);a,/3 > 0}, T3 = {iffi(/3a;);a,/3 > 0} and Ti = {i/;(/3x);a,/3 > 0}. 

(2) It is clear that if / is an LWMY function, then the functions / ^ and x n> 
i/(/3a;), a, /3 > 0, are LWMY functions. 

(3) The image of J-2 by the map / n- /^^ is T^. The functions x 1— > a/x and fs are 
involutive. 

In the sequel, we focus on the three new cases: cither / = /i, f = gi or f ~ fl and in 
each case, we determine the laws of the related random variables. 
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2.1. The cases f = gi and f = fi 

(a) Recall the definitions of the gamma distribution 7(A,c)(da;) = 
jYXT2'''^~^e~^^l(o,oo)(2;) da;(A,c > 0) and the beta distribution Beta(a,&)(da;) = 

p/'^p.A x°~^(l — x)^~^l{o<a:<i}dx(Q,b > 0). Considcr (sec, e.g., [13], or [14] and the 
references therein) the Kummer distribution of type 2: 

K'-^Ha, h, c) := a{a, h, c)x''~^{l + a;)-"-''e-=^l(o,oo) {x) Ax, a, c> 0, 6 £ M, (2.14) 

where a{a,h,c) is a normalizing constant. 

Associated with a couple {X, Y) of positive random variables, consider 

(t/,F):=r,,(X,y)=(-^^,-^--^^). (2.15) 

In Theorems 2.4 and 2.6 below, we suppose that all random variables have positive 
and twice differentiable densities. 

First, we consider the case f — fi- We determine the distributions of X and Y such 
that /i is an LWMY function associated with {X, Y) . 

Theorem 2.4. (1) Consider two positive and independent random variables X and Y . 
The random variables U and V defined by (2.15) are independent if and only if the 
densities of Y and X are, respectively, 

P'-^y^ = fW§)^^ e-^)^-ie-^l{,>o}, (2.16) 

px(x) = a(a + b, c, -a)e-(''+'')^(l - c-^y"-^ 
xexp(^-c ^°^_^ jl{^>o}, 

where a, b and c are constants such that a,b,c>0 and a(a + b,c, —a) is the constant 
from equation (2.14). Thus, the law of Y is the image of the Beta(a,6) distribution by 
the transformation z £ (0,1) H- — logz G (0,oo), while the law of the variable fi{X) is 
K^'^\a + b,—b,c) (see equation (2.14)). 

(2) // (1) holds, thenU^K'^^\a,b,c) and V ^-f{b,c). 

The proof of Theorem 2.4 will be given in Section 4. 

Remark 2.5. Since gi = fi^, Remark 2.3 and Theorem 2.4 imply that the random 
variables associated with the LWMY function gi are the random variables U and V 
distributed as in item 2 of Theorem 2.4. 

(b) As suggested by identities (2.16) and (2.17), it is possible to simplify the statement 
of Theorem 2.4. Since Tg-^ = T7^, we have 

(X.,..„,.,.,.(,o.(l±^),,o.(l±^)-,„S.(i±^)). ...S, 
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As (2.18) shows, it is useful to introduce 

Obviously, the correspondence {U, V) i— > (C/', V) is one-to-one: 

(u V) ( ^'^^ v'iv' + m-u') \ 

Furthermore, {X,Y) can be easily expressed in terms of {U',V'): 

X = log(l + l/F') and Y = -\ogU'. (2.21) 

Since it is easy to determine the density function of 4>{^) knowing the density function of 
a random variable ^, where (j) is differentiable and bijective. Theorem 2.4 and its analog 
related to / = gi (see Remark 2.5) are equivalent to Theorem 2.6 below. 

Theorem 2.6. (a) Let U' and V' be two positive and independent random variables. 
The random variables U and V defined by (2.20) are independent if only if there exist 
some constants a, b, c such that 

[/' - Beta(a, &) and V' ^ K^^\a + b,-b,c). (2.22) 

If one of these equivalent conditions holds, then U ^ K^ ''{a^b^p) and V ~7(6, c). 

(b) Let U and V be two positive and independent random variables. The random vari- 
ables U' and V defined by (2.19) are independent if only if there exist some constants a, 
b, c .such that 

U^K'''^\a,b,c) and y- 7(6,0). (2.23) 

Under (2.23), C/' ~ Beta(a, fe) and V ^ K'^^^a + b,-b,c). 

We now formulate a simple consequence of Theorem 2.6. 

Theorem 2.7. For any a,b,c> 0, the transformation {u,v) i— > ( — i+Wtt" ,u + v) maps 

the probability measure K^'^'{a,b,c) ® ^{b,c) to the probability measure Beta(a,&)(8> 
K^'^'{a + b,—b,c). In particular, 

K^^'l {a,b,c) * j{b,c) = K^^\a + b, -b,c). (2.24) 

Remark 2.8. Note that (2.24) may be regarded as an analog of (1.2). 

2.2. The case / = /; 

Recall that /^ has been defined by (2.13). Due to the form oi fg , a change of variables 
allows us to simplify the search for independent random variables X and Y such that 
the two components of Tf {X, Y) are independent. 
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For any decreasing and bijective function / : (0, cxd) — > (0, oo), we define 

7(a;)=cxp{-/(-loga;)}, xe(0,f), (2.25) 

Tpix,y)^(^f{xy),jl^^, x,ye(0,f). (2.26) 

Observe that / is one-to-one and onto from (0, 1) to (0, 1), T^ is one-to-one and onto 
from (0, 1)2 to (0, 1)^ and 

(Tf)-i=T;^i. (2.27) 

Definition 2.9. Let X and Y be two independent and (0,l)-valued random variables. 
We say that a decreasing and bijective function / : (0, 1) — > (0, 1) is a multiplicative LWMY 
function with respect to {X,Y) if the random variables [/'" := f{XY) and V™" := /(xy) 
are independent. 

Remark 2.10. For any random vector {X,Y) in (0,00)^, we consider X' = e~^ and 
Y' = e^^. Then, / is an LWMY function with respect to {X,Y) if and only if / is 
a multipUcative LWMY function with respect to {X',Y'). 



The change of variable x' = e ^ is very convenient since the function 

l-x 

+ (<5-l)x' 



cj,s{x):=fi{x) = -^-^^—^, xe{OA) (2.28) 



is homographic. 

Note that /| : (0, 1) -^ (0, 1) is bijective, decreasing and equal to its inverse. First, let 
us determine the distribution of the couple (X',Y') of random variables such that (ps is 
a multiplicative LWMY function with respect to {X',Y'). 

For a, 6, a > and c G M, consider the probability measure 

/3„(a, 6; c){dx) = ka{a, b; c)x''~\l - x)''-\ax + 1 - a;)^l(o,i)(x-) dx. (2.29) 

Note that if c = 0, then I3a{a, b; c) = Beta(a, b). 

Theorem 2.11. Let X' and Y' be two independent random variables valued in (0,1). 
Consider 



l + {6-l)X'Y'' l + {6-l)X' l-X'Y' 

for fixed 5 > Q. 

Then, C/™ and V^™ are independent if and only if there exist a,b,X> such that 

X' r^Ps{a + b,X;-X-b), y' - Beta(a,6). (2.30) 
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If this condition holds, then 

U"'r^PsiX + b,a;-a-b), F™ ~ Bcta(A,6). (2.31) 

In the case (5 = 1, Theorem 2.11 takes a very shuple form. 

Proposition 2.12. Let X' and Y' he two independent random variables valued in (0, 1). 
Then, 

U"^ = l-X'Y', V"'= ^^Sx'Y, 

are independent if and only if there exist a, 6, A > such that 

X'-Bcta(a + 6,A) and F' - Bcta(a,6). 
// one of these conditions holds, then [/™ ~ Bcta(A + b, a) and V^^ ~ Bcta( A, b) . 

Remark 2.13. When X' -- Bcta(a + 6, A) and Y' ~Bcta(a,6), it can be proven that U™ 
and T^™ are independent using the well-known property that if Z and Z' are independent 
with Z ~ 7(0, 1) and Z' ^ 7(6, 1), then R := „jf„, and Z + Z' are independent with 
R - Beta(a, b) and Z + Z' '^j{a + b, 1) (see, e.g., [17]). 

According to Remark 2.10, fg is an LWMY function with respect to {X, Y) if and only 
if (j)s is a multiplicative LWMY function with respect to {X' , Y') = {e^-^ , e~^) . Therefore, 
a classical change of variables allows us to deduce that Theorem 2.11 is equivalent to 
Theorem 2.14 below. 

Theorem 2.14. (1) Consider two positive and independent random variables X and Y. 
The random variables U = fg{X + Y), V = fgiX) — fg{X + Y) are independent if and 
only if the densities of Y and X are, respectively, 

P'-^y^ = fW§}^^ c-^)^-^c-^l{.>0}, (2.32) 

-~x\ — X — b 

(2.33) 



px{x) = ks{a + b, A, -A - 6)e-(°+'')^((5e~^ + 1 - e 
x(l-e--)^-il,>o. 



where a,b> 0, A G M and ks{a + b, X,~X~ b) is the normalizing factor (see (2.29)). Thus, 
e~^ is Beta(a, b)- distributed and e^^ is (3s{a + 6, A, —A — b)- distributed. 
(2) // (1) holds, then the densities of U and V are, respectively, 

pu{u) = ks{X + b,a; -a - 6)e-"(^+'''(l - c^")"-^ 

(2.34) 

x(l + (,5-l)e-")-'^-''l„>o, 
Pv{v)^c-^-{l-c-'')''-H,,yo. (2.35) 

We omit the proof of Theorem 2.14 since it is similar to that of Theorem 2.4. 
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3. The set of all possible "smooth" LWMY functions 

The following theorem gives a functional equation linking LWMY functions to the related 
densities. 

Theorem 3.1. Let X and Y be two independent and positive random variables whose 
densities px and py are positive and twice differentiable. Define (px = ^ogpx and 
(f>Y = logpy. Consider a decreasing function /:(0,cxd) i— >• (0,oo), three times differen- 
tiable. Then, f is a LWMY function with respect to {X,Y) if and only if 

J"{x) ,1 , ( ^ 1 



fix) ^''"'' ^ 'yf'{x) f'ix + y) 

rM , 2{f"{x)r-f "{x)f{x) 

fix) fix. 



- VvW-fTT:^ + TTTAi =^' x,y>Q. (3.1) 



Proof. Let g = f'^ and {U,V) = Tf{X,Y). By formula (2.9), {X,Y) = Tg{U,V). X 
and Y being independent, the density of {U, V) is 

V(u,v){u,v) = px{g{u + v))pY{g{u) ~ g{u + v))\J{u,v)\lu,v>o, (3.2) 

where J is the Jacobian of the transformation Tf. We get \J(u,v)\ ~ g'{u + v)g'{u), and 
then 

P{u,v){iJ; v) = px{g{u + v))pY{g{u) - g{u + v))g'{u + v)g'{u). (3.3) 

The variables U and V are independent if and only if the function H — \ogP[u,v) 
satisfies q^q^ = 0. By equation (3.3) we obtain 

1^ = cl>';,{x)[g'{f{x))Y + M^)g'\f{x)) 

-ct^'i-{y)g'{f{x))[g'{f{x + y))-g'{f{x))] 
9'" 9' - id"? 

{g'Y 

where x = g{u + v) and y = g{u) — g{u + v). Differentiating three times the relation 
gifix)) = X, we obtain g"if{x)) - -^ and g"'{f{x)) - - riMH^^t^^n^ _ As a re- 
suit, 

9"'9'-{9")' ,f ,, 2f"{xf-f"'{x)f'{x) 

Therefore, ^fj = leads to (3.1). D 

We restrict ourselves to smooth LWMY functions /, that is, those satisfying 

/ : (0, oo) —5- (0, oo) is bijective and decreasing, (3.6) 



4.'y(y)g"{f{x)) + ^ ' ,;/ ' (fix)), (3.4) 
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f is three times differentiable, 
F{x) = Y^ a„a;" Vx > 0, 



(3.7) 
(3.8) 



?i>i 



where F := 1//'. 

According to (3.6), /'(0+) = — oo. This implies that -F(0+) = and explains why the 
series in (3.8) starts with n~l. 

The goal of this section is to prove half of Theorem 2.2: if / is a smooth LWMY 
function, then / belongs to one of the four classes J-i, . . . , J-4 introduced in Remark 2.3. 
First, in Theorem 3.2, we characterize all possible functions F. Second, we determine the 
associated functions / (see Theorem 3.9). 

Theorem 3.2. Suppose that f is a smooth LWMY function and the assumptions of 
Theorem 3.1 are satisfied. 



1. IfF'iO^ 

F{x) 

2. IfF'(0+)^0, then 



0, then 02 < and 

coshl X 



a\ 



12a4 

02 



a^x 



if aA < 0, 
otherwise. 



(3.9) 



Fix) 



0102 



coshl X 



3a3 
L aix + a2X 




Ol 

603 

2 



603 
ai 

sinhl X 



if aiaz >0, 
otherwise. 



(3.10) 



that 



Remark 3.3. Unsurprisingly, the case F{x) = a^x^ corresponds to f[x) = — - 
is, the case considered by Matsumoto and Yor, and Letac and Wesolowski. Thus, under 
stronger assumptions, we retrieve the result of Letac and Wesolowski. Indeed, writing 
the functional equation of Theorem 3.1 with / : a; i— > \/x gives 

0x(^) + -<t^'x{^) + My)\{^^ - i^ + vf) -My) + 4 = 0. 

X X'^ X x^ 

We then solve this differential equation and find that the laws of X and Y are necessarily 
GIG and gamma, respectively. We omit the details. 

Throughout this subsection, we suppose that / satisfies (3.6)-(3.8) and that the as- 
sumptions of Theorem 3.1 are fulfilled. To simplify the statement of results below, we do 
not repeat these conditions. 

Recall that 0y is the logarithm of the density of Y . Let us introduce 



h:=(j)Y- 



(3.11) 
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Lemma 3.4. 

1. There exists a function A : (0, cxd) — > M such that 

, \(x) - h(y)F'(x) , , , , 

Fix + y)= ' ' ^,^^J^^ ^ ' +Fix). (3.12) 



2. F satisfies 



\{0+)-h{y)F'(Q+) 



Remark 3.5. Suppose that we have been able to determine F. Then, h ~ (p'y solves 
the linear ordinary differential equation (3.13) and can therefore be determined. The 
remaining function (j}x is obtained by solving equation (3.1). 

Proof of Lemma 3.4. Using (3.11) and F = 1//' in equation (3.1), we obtain 

c{x) = Hy)^ + h'iy)-^iFix + y) - Fix)), 
b (x) b (x) 

where c{x) depends only on x. Multiplying both sides by F{x) and taking the y-derivative 
leads to 

= F'{x)h'{y) + [F{x + y)- F{x)]h"{y) + h'{y)F'{x + y). 

Fix X > 0. Then, 6{y) :— F{x + y) is a solution of the differential equation in y 

= F'{x)h'{y) + {9{y) - F{x))h"iy) + h'iy)9'{y). (3.14) 

A solution of the related homogeneous equation in y is /pf-r, where p is a constant. It is 
easy to prove that y i~-> —F'{x)h{y) + F{x)h'{y) solves (3.14). Thus, the general solution 
of (3.14) is 

0{y) = )^[^(^) - F\x)h{y)+F{x)h'{y)]. 

Since e{y)=F{x + y), (3.12) follows. 

According to (3.8), F{0+) and F'{0^) exist. Therefore, taking the limit a; — > 0+ 
in (3.12) implies both the existence of A(0+) and relation (3.13). D 

The following lemma shows that the function F (and thus /) solves a self-contained 
equation in which h, and thereby the densities of X and Y , are not involved. 

Lemma 3.6. F solves the delay equations 
X{0+)-h{y)F'{0+) 
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F'(^^ + y)^EM±pMlF{a: + y)-F{x)]-F'{x) {x,y>0). (3.16) 

Proof. By (3.13), wc have 

A(0+) - h{y)F'{0+) 



h'{y) 



F{y) 



Equation (3.15) then foUows by rewriting equation (3.12) and replaeing h'(y) with the 
expression above. 

We differentiate (3.15) in y and use the fact that A(0+) - h{y)F'{0+) == h'{y)F{y) to 
obtain 

By (3.12), we have ""^.7^,'^'" = "'"47"'" and this give, (3.16). O 

Remark 3.7. We can see (3.16) as a scalar neutral delay differential equation. Indeed, 
set t — X + y and consider j/ > as a fixed parameter. Then, (3.16) becomes 

F'[t) = a{F{t)~F{t-y))-F\t-y), t>y, (3.17) 

where a := ^''^y^Jj'"+^ Replacing F{t) in (3.17) with e"*G(t) leads to 

G'(i)+e-"«G'(t-y)+2ae-''^G(i-zj) = 0, t>y. (3.18) 

Equation (3.18) is called a neutral delay differential equation (see. e.g.. Section 6.1, in [6]). 
These equations have been intensively studied, but the authors have only focused on the 
asymptotic behavior of the solution as t ^- oo. Unfortunately, these results do not help 
to solve explicitly either (3.16) or (3.18). 

Lemma 3.8. For all integers fc > and 1>1, we have 

i-i 

^ (/ - 2m + l)Gfl„_|_i_|_j.a;_,„+i+fea„ = (l - 2)(fc + l)ak+iai + aia;+feG;+j.,(3.19) 

Gfc+gOfc+sai = (fc + l)afc+ia3, (3.20) 

2Gf+4afc+4ai + Gf+3afe+3a2 - C^+20'k+2a3 - 2(fc + l)ak+ia4 = 0, (3.21) 

where Cf, ^ , "' , , . 

Proof. Obviously, the equation (3.16) is equivalent to 

F'{x + y)F{y) - F'{y)F{x + y)- F\y)F{x) 

(3 22) 
- F{y)F\x) + F'{Q+)F{x + y)~ F'{0+)F{x). 
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Using the asymptotic expansion (3.8) of F, we can develop each term in (3.22) as a series 
with respect to a; and y. Then, identifying the series on the right-hand side and the 
left-hand side, we get (3.19)-(3.21). The details are provided in the Appendix. D 



Proof of Theorem 3.2. We will only prove item 1; the proof of item 2 is similar. 

Since ai = F'(O-i-) ~ 0, we necessarily have 02 ^ 0. Indeed, if 02 = 0, then, by (3.21) 
with fc = 1, we would have —3a| — 40204 = 0, that is, 03 = 0. Again using (3.21) with A; = 3 
would imply that 04 = and finally that Ofc = for every A: > 0, which is a contradiction 
because, by definition, F ~ 1//' does not vanish. 

So, we have oi = and 02 ^ 0. Equation (3.20) with fc = 1 reads 40401 = 20203, which 
implies that 03 ~ 0. Applying (3.20) to fc = 2n provides, by induction on n, 02n+i = for 
every n > 0. 

Therefore, equation (3.21) reduces to (fc4-3)(fc-|~2)(fc+ 1)0^+302 = 12(fc+ 1)0^+104, fc > 



0, that is, Ofe+3 = ^ (k+3)(k+2) "'k+i- This leads to 



a-ik 



fUoiV ^ 2 



Ur 



(2fc-)! 



a2, 



fc>l. 



(3.23) 



Then, F{x) =022;^ if 04 = 0, and if 04 7^ 0, we have 

fc-i , 



fc>i ^ 



/1204 



02 



:a2X 



2 k 



If O4O2 < 0, then F{x) ~ 
is impossible since F(x) 



-12a4 



7r^fcos(a; 
= l/f'{x) < 0. Consequently, 



(2fc)! 
1]. This implies F{2n 



-12a4 



0, which 



F{x) = -p 
D04 



cosh I X 



1204 \ _^ 



02 



D 



Now, in each case of Theorem 3.2, we compute the function / associated with F via the 
relation F = 1//'. Wc do not detail the calculations since they reduce to getting a good 
primitive of l/F. Recall that we restrict ourselves to functions / satisfying (3.6)-(3.8) 
and work under the assumptions of Theorem 3.1. 



Theorem 3.9. 



1. IfF{x)^a2X^, then fix) ^^. 



2. //F(a;)=a(cosh/3.T-l),a,/3>0, then f{x) ^ ^f^i^x). 

3. IfFix)^aix + a2X^, then fix) ^ -^gii^^x). 

4. // 



Fix)- 



0102 
303 





sinhl X 
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then 

ff \ 1 , fe^^ + S-1 

^here a= ^, /3 ^ ^ and j ^ aiy^. 

4. Proof of Theorem 2.4 

Recall that (py = logpy, h — (f)Y and F'(O-i-) = 0. It is easy to deduce from (3.13) that 
there exist constants A and ci such that h{y) ~ ^f{y) + ci , that is, h{y) — —;^z^ + ci — A. 
This implies the existence of a constant d such that (f)Y (y) — A log(e*' — 1) + (ci — X)y + d. 
Setting M = e'*, we have, by integration, for all y > 0, 

PY{y)^M{l-c-yfc''^y. (4.1) 

To give more information on the normalizing constant M , we observe, for a ~ — ci and 
6 = A + 1, that 

/■oo /•! 

f M{\ - e-yf-^e-'^y dy = M / (1 - uf-^"-^ du, 
Jo 

which implies that a > 0, 6 > and M = Y-uycib) ■ "^^^^ proves (2.16). 

To find the density of X, we return to equation (3.1) and compute each of its terms. 
We have f{x) = j^, f"{x) = ^^ and /'"(x) = -2!l±i£fl+el ^^ that ^^ = 

"(Ir^-i)^^^ ^'^'^ ^IW) ^ ~S^- Calculations yield 



2(/"(:^))'-/"'(a;)/'(x) 



/'(x)2 (e^~l)2' 

Moreover, 

-,2a; 



(4.2) 



,, , J"{^) , .,, J rW A (ci~A)e2--ci 

~ ^-^^^7R + ^^^"^ \J¥Ty) " V = (0-1)2 • (4-3) 

Equation (3.1) can then be written, using (4.2) and (4.3), 

w,"^ ^^«" + l^' r ^ (ci-A-l)e2— ci-l 

Then, ft-o '■— 4>'x solves 

^// N e^ + 1, / ^ (ci-A-l)e2^-ci-l .^ ^. 

'^o(^) + ^^3T^o(^) = ^ J^^2 ■ (4-4) 
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Note that x i— > . . Jf, , . solves (4.4) with the right-hand side equal to 0, and x H> 
^'^^~ 4sinii4'3;'/2T — ^^ ^ particular solution of (4.4). Therefore, the solution of (4.4) is 



h{x) 



(ci-A-l)e^ + (ci + l)c-^ + A' 



4sinh^(a;/2) 
for some constant K . This implies that 

(2ci-A + A')e^ (A + 2)e=^ 



6^(x)=ci + l + 



(e^-l)2 



1 



As a consequence, there exists a constant S such that 



cl>x{x) = (ci + l)x - ^^""^ ,'^V^'''^ - (A + 2) log(e^ -1) + S. 



Thus, pxix) = iVe(=i+i)^(e^ - l)-^-2exp(-2£^^^)l{^>o}. Recah that a = -ci and 
6 = A + 1. With c = 2ci — A + A', we get (2.17). More information on the constant N is 
obtained by observing that if we set V' = fi{X) = -pr^, then the density of V' is 



fy,(w) = Niw + l)-'^u.'^+''-i exp{-c«;}l{„>o}, 

that is, the law of V is Ar(^)(a + 6, — 6,c) (see equation (2.14)). 

We have g'i{u) — — ^tu+i) - ^ computation of a Jacobian, together with (2.16) 
and (2.17), implies, for u,v>0, that 



P{uy){u,v) =pxilog 



U + V + 1 



u + v 



PY log 



{u+l){u + v) 

u{u + f + 1) 



u{u + l){u + v){u + tl + 1) 



We then get that p([/y)(u,w) is the product of a function of u and a function of w, and 
this gives item 2 of Theorem 2.4. 



Appendix 



Proof of Lemma 3.8. We have 



F'{x + y)F{y)^^x'' ^ na„a™C;^_iy 

fc>0 m>0,n>l+fc 



k n+m—l — k 



Setting l~m + n — \ — k for fixed m gives 



F'{x + y)F{y) 



: ^ x*^?/' ^ (/ - m + 1 + fc)Cf_,„+j,a,_,„+i+fca™. (5.5) 

A;>0,/>0 m=0 
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By the same method, we have 

F'{y)F{x + y)= ^ x%^ ij^ mCtm+k+iai-m+i+kam]. (5.6) 

k>Od>0 \m=0 / 

As for the two other terms of (3.22), we get 

F'{y)Fix)= Y. akai+i{l + l)x''y\ (5.7) 

fc>0,/>0 

F'{x)Fiy)= Y, ak+iai{k + l)x''yK (5.8) 

k>0,l>0 

Consequently, 

F'{0+)F{x + y) = ai 5Z ""(^^ + vT = «i 5Z (^i+kCt+k^^^y', (5.9) 

Tl>0 fc,i>0 

F'{Q+)F{x) = ai ^ akx^. (5.10) 

A:>0 

Identifying the coefficient of x^y^ in (3.22) and using (5.5)-(5.10), we have, for fc > 
and I > 0, 

I 
^ (/ - 7Ti + 1 + fc)Cf_„j_|_i,a;_„i+i+fca„j = -(/ + l)akai+i - (fc + l)ak+iai 

l+l 

+ Y mC^_^.^^f,+iai-«i+i+kam (5.11) 

+ aiaij^kCi^j. — aittfel/^o. 

Note that if Z = 0, then both sides of (5.11) vanish. Therefore, we may suppose in the 
sequel that / > 1 . 

For m = / + 1, we have mC/L,„+fe+ia(-m+i+fcam = (^ + ^)o-kO-i+i- Thus, equation (5.11) 
reads 

I 
^ (Z - TO + 1 + fc)C/l,„+j,a,_„+i+fca„ 

I 

= -(fc + l)afe+ia; + ^ TOCf_„j+;,^ia/_,„+i+fca„ (5.12) 



m=0 



+ aia/+feC,^fe. 
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However, via a calculation involving the definition, we find that 

(/ - m + 1 + k)Ctm+k - "^Q-m+i+fc = (^ - 2m + l)Ct„.+i+k. 

so equation (5.12) is equivalent to 

I 
^ (? - 2m + l)C/LTO+i+fca;_m+i+fcam = -(fc + l)afc+iai + aiai+kCf^^- (5.13) 

For m = l, we have (/ - 2to+ l)C/l^^^_|_j,a;_m+i+fea,„ = (1 - 0(fc + l)afc+ia;- Conse- 
quently, equation (5.13) may be written as 

i-i 

22 (^-2TO + l)C/l,„+i+fca;-m+i+fcam-(^-l)(fc + l)afc+ia; = -{k + l)ak+iai+aiai+kCl^k, 

which implies (3.19). 

(3.20) and (3.21) follow by applying (3.19) to / = 3 and / = 4, respectively. D 
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